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Abstract 

The persistent currents driven by the pure Aharonov-Bohm type 
magnetic field in mesoscopic normal metal or semiconducting cylinders 
are studied. A two-dimensional (2D) Fermi surfaces are characterized 
by four parameters. Several conditions for the coherence and enhance- 
ment of currents are discussed. These results are then generalized to 
a three-dimensional (3D) thin-walled cylinder to show that under cer- 
tain geometric conditions on the Fermi surface, a novel effect - the 
appearance of spontaneous currents is predicted. 
PACS numbers: 71.30.+h, 72.10.-d 

Keywords: mesoscopic cylinder, Fermi surface, spontaneous currents. 



1 Introduction 

Quantum coherence in mesoscopic systems arose a great interest in the last 
years because it leads to many interesting phenomena such as persistent 
currents, conductance oscillations etc. 0. In a series of papers |2|, || we 
investigated persistent currents in a collection of isolated mesoscopic rings 
stacked along one axis. The current I(<p) induced by the magnetic flux <j), due 
to the Aharonov-Bohm effect, was therein calculated under the assumption, 
that the electrons from different rings interact via the magnetostatic (current- 
current) interaction and the interaction was considered in the mean field 
approximation. It was shown that due to the selfconsistent equation for 
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the current, in some cases, a novel effect - the spontaneous selfsustaining 
currents at the zero external flux - can be predicted. A similar effect, due to 
the Aharonov-Casher effect was also recently proposed Q]. 

In the present paper we investigate persistent currents in 2D and 3D cylin- 
ders made of a clean normal metal or a semiconductor (ballistic regime). It 
is known that the magnitude of the currents in 2D systems depends strongly 
on the correlations of currents from different channels f|. To study this 
correlation we calculate the currents in systems with different shapes of the 
Fermi surface (FS). We also search for the possibility of obtaining sponta- 
neous currents for the most favorable FS shapes. 

2 Persistent Currents in a Single Channel 

Let us consider a 2D normal metal or semiconductor system of a cylindrical 
geometry. The circumference and height of the cylinder are denoted by L x 
and L y respectively. We neglect the width of its wall, assuming it is of the 
order of - the Fermi wavelength. We also assume that the magnetic flux 
threads the cylinder axially and is confined to its center, so that the electrons 
move in a field-free space. Here we use the cylindrical coordinate system, in 
which the height coordinate is y and the azimuthal angle is G G [0; 2ir). It 
is convenient to replace G by x = so that x varies between and L x . 

For the wave function of an electron in the cylinder we apply periodic 
boundary conditions in the azimuthal direction and infinite potential barrier 
along the y axis. In presence of the magnetic field <p the wave function obeys 
the following boundary conditions || 



where (po = ~ denotes the flux quantum. 

Solving the Schrodinger equation in the free electron approximation we 
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get quantized energy levels 



M*) = ^[*2M) + W1> ( 2 ) 

where k x (q, <fi) and k y (r) are components of the wave vector k in the azimuthal 
and y direction, respectively, 

k x (q,<P) = ^(q + ^), g = 0,±l,±2,..., (3) 

TTX 

ky(r) = —, r = 0,l,2,.... (4) 

y 

In our 2D system we are interested in the net current along the azimuthal 
direction. The Fermi surface (FS) is in this case 2D. Apart from momentum 
in the k x direction the electrons can move in the y direction. The collection 
of them, having the same k y or r quantum number will be called a channel. 
Each channel is a semi one-dimensional (ID) system with a characteristic 
radius K x , which is the maximal possible value of the x component of the 
wave vector \k x \ < K x . This radius as well as the number of channels M 
depend on the shape of the 2D FS. 

The current carried by a single electron of the energy E a is given by 

W) = — (5) 

where a is an arbitrary set of indices. In case of the energy (0) it is a 
linear function of 0. The total current in the 2D system is the sum over 
the contributions I a of all states, where the occupation probability is the 
Fermi-Dirac (FD) function. 

In the present and in the the next chapter we assume that the temperature 
is zero. We first recall the current I r {4>) corresponding to a single channel r 
with N electrons. It is a sum over q of currents I q r{4>) corresponding to the 
first N lowest lying states E qr (<f). In the following we use the Fourier series 
expansion formula [EL [7| for this current 

W) = y2~r cos ( lK ^ L ^ sin \~~r~ J ' ( 6 ) 

1=1 7r/ V 0o J 
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where 



evp ehK x . . 



is the amplitude of the current in the ID system and K x = ^ is the ID 
radius of the FS corresponding to the channel r {vp is the velocity of electrons 
at the FS). The total current for the single channel r has a sawtooth shape 
with the period 1 in ^ and the amplitude inversely proportional to the 
circumference of the ring. The current for a system with an even number of 
electrons has a discontinuous jump up at <fi = whereas in the case of an odd 
number of electrons at 4^ = 0.5. More generally if the number of electrons 
in the system varies between even and odd as a function of (j), the current 
assumes one of the above values corresponding to the number of electrons at 
the particular value of the flux |5], §]. In the next chapter we show that this 
property to some extent "survives" the transition to a 2D system. 



3 The Currents Depending on the Shape of 
the Fermi Surface 

The total current in the 2D case depends on the strength of the interchannel 
correlations. It has been already shown |J that if the FS is spherical then, for 
short cylinders L x > L y the channel currents add without phase correlation, 
whereas for long cylinders there exist some phase correlation and the total 
current is bigger. A large phase correlation among channel currents means 
that the increase of the flux <f> results in an almost simultaneous cross of the 
FS by the large number of channels. The most favorable situation takes place 
if the separation between the last occupied level and the FS from channel to 
channel is nearly the same. There exists then a perfect correlation among 
the channel currents because the M levels cross the FS simultaneously while 
the flux is changed by one fluxoid - we get then the largest amplitude of the 
total current. 

The phase correlation and the value of the total current depend on the 
shape of the FS. In case of the free electron approximation considered in 
chapter 2 the FS is circular 
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K 2 F = K> + K> r (8) 

It is characterized by a single parameter Kp, specifying its volume. The 
radius K x of the FS in the given channel r can be easily calculated from @. 

However for the interacting electrons moving in the crystal lattice the 
dispersion relation E(k) is in general a complicated function of the wave 
vector k, leading to different shapes of the FS. For example, in the tight 
binding approximation, for 3D systems we can get, depending on the crystal 
symmetry and on the filling factor, the FS being the sphere, the cube, the 
cube with rounded corners, octahedron etc. 15). The FS can be in addition 
differently oriented with respect to k a (a = x,y,z) axes in the reciprocal 
lattice. For 2D systems or for 3D systems with 2D conduction (an example 
of such systems are high T c superconductors in a normal state) the FS can 
change from a circle to a square [ID, IT] via different other shapes. 



In the following we study the influence of the shape of the FS on persistent 
currents in multichannel systems changing its shape in the following way 



K x 



The shape of this Fermi surface depends on four positive parameters: Kp, n, 
a and f3. Generally speaking the parameter n controls the convexity of the 
FS, whereas a and j3 measure its curvature at the places where it crosses k x 
or k y axes. To understand how the above parameters influence FS and the 
corresponding current in this chapter we consider three cases in which some 
of the parameters are kept constant while the other are varied. 

(i) Kp = const., n = 2 and a, j3 are varied: 
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K 2 K 

In this case the FS is elliptical with its major axis along k x or k y 
depending on whether a > (3 or (3 > a, respectively. The curvature of 
the FS at the point, where it crosses the k x axis is equal to ^K F . We 
assume in addition that aj3 = 1 because it ensures the independence 
of the volume under the FS on a and j3. Under this assumption the 
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number of states under the FS maintains independent on a (and (3) in 
the limit of Kp — > oo whereas for finite Kp it depends only slightly 
due to the effects on the edge of the FS. 

Replacing K y by its actual value at the r-th channel (^) the Eq. (PH|) 
yields for the ID radius K x of the FS in this channel 



K x = aK F 



\ 



1 - 



'kyijry 
BK F 



(11) 



The current in the whole 2D system is the sum of M + 1 single channel 
currents @ taken with the appropriate K x 



M 



r=0l=0 ' U 



\ 




COS 



lKpL x a 



\ 



1 - 




ehK F 
mL x 



sm 



;i2) 



'2-Klf 



where I 

For the circular Fermi surface (a = (3 = 1) in a 44-channel system for 

□ h 2 K 2 
L x = 1000A, L y = 100A and assuming E F = „ F = 7eV (typical 



for Cu) the characteristic is shown in Fig.l (solid line). It is a insert Fig.l 
ragged function, built of intervals of the linear decrease and discon- 
tinuous jumps up of the current. Each discontinuous increase of the 
current value in Fig.l corresponds to a single electron leaving or com- 
ing under the Fermi surface. This structure of the current function is 
characteristic for T = and remains the same in other cases in this 
chapter. 

If the parameter a increases then the circular FS turns into an ellipse 
with the major axis along the k x -&xis. An increase of a up to 10 (Fig.l - 
dashed line) results in the decrease of the number of channels (M ~ (3) 
and the current value. In this case we don't observe any correlation 
between currents from different channels. An increase of the parameter 
(3 is equivalent to turning the major axis of the ellipse along /c y -axis. If 
(3 increases up to 10, then a substantial enhancement of the current is 
observed (Fig.l - dotted line). This is because the perpendicular to the 
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k x axis part of the FS is very flat (its radius of curvature is proportional 
to causing an enhancement of the current correlation. 

It is worth noticing that the current correlation change similarly for 
the circular FS (a — f3 — 1) when the circumference and the height 
of the cylinder are changed in such a way that L x L y = const.. In this 
case an increase of the cylinder height produces the same enhancement 
of current correlation as the increase of (3. The total current however 
is even more enhanced because, at the same time, the circumference of 
the cylinder decreases causing the increase of I (cf (|i~2D). 

(ii) K F = const., a — (3 = 1, n is varied: 



K F = \K x \ n + K™. (13) 

Here we consider n to be an arbitrary positive real number. For n = 1 

the Fermi surface is of a triangular shape (Fig.2). With increasing n it insert Fig. 2 

becomes convex and changes from triangular through circular (n = 2) 

to rectangular for n — > oo. On the other hand, decreasing n below 

n — 1 yields concave Fermi surfaces shown in Fig. 3 . Such FS are insert Fig. 3 

frequently observed in HTSC materials. In our numerical calculations 

we adjust Kp to n in such a way that the volume under the FS remains 

constant. 

In the present case K x is 



k~ = Kp r 



'k y (r)\ n 



Inserting this into (|5[) and (|7]) yields 



(14) 



M oo 



M0 = ££?t 



r=0 1=0 



Til \ 



1 - 



' kyir) 
Kp 



cos 



lK F L a 



\ 



1 - 



Kir) ' 
Kp 



sin 



'2nlf 



(15) 

The case with n = 0.4 corresponds to the FS from Fig. 3. It is the 
case where the number of channels is greater compared to the circular 
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FS (n = 2). Taking again L x = WOO A, L y = 100 A and E F = 7eV 

the I— 4> characteristic for n = 0.4 is shown in Fig.4 (dotted line), insert Fig. 4 

We observe here a moderate increase of the current value compared 

to the circular FS (solid line). For n = 1 the FS is triangular and 

in the I—<p characteristic (Fig.4 - dashed line) there is a substantial 

increase of correlation of phases from different channels and the greater 

current value. It is because all the states at the FS leave (or get in) 

the FS simultaneously. More generally, in case of n = 1, there is a 

full correlation and enhancement of the current provided the following 

geometrical condition 



^aL x = fx(3L y , (16) 

where \jl is an arbitrary positive integer, is obeyed. In case of the dashed 
line of Fig.4, fi = 5. 



However, if the relation (If:) do not hold the 1(0) current can be much 

smaller. Therefore the geometrical amplification of the current in case 

of n = 1 is a matter of a careful adjustment of the cylinder dimensions 

L x and L y in such a way that, for a given a and (3, the condition (pf) 

is fulfilled. Such a dependence on a cylinder dimensions does not take 

place if we consider persistent currents for the systems with n > 2. For 

n = 5 we have a case in between circular and rectangular FS. In the I-cp 

characteristic (Fig. 5 dotted line) some correlation among the channel insert Fig. 5 

currents exists what results in an increase of the current amplitude 

compared to the circular FS (solid line). The case with n —>■ oo, i.e. 

rectangular FS, is the most suitable to obtain a big current amplitude. 

The I—(f) characteristic is shown in Fig.5 (dashed line). We observe 

here a perfect coherence (all the terminal states simultaneously leave 

or get into the Fermi surface) and great increase of the current value. 

The amplitude of the current is equal to the number of channels M in 

Jo-units. Thus in general, with the departure of n from 2, we observe 

the increasing amplitude of the M-channel current. 

(iii) a f3, n = 2, K F is varied: 

In the previous two cases we have changed the shape of the FS leav- 
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ing the area under it invariant. As a result of that some currents 
were enhanced and some were suppressed compared to the spherical 
FS. Changing parameters a, (3 and n we could not, however, predict 
whether the current is positive or negative in a given interval. This im- 
portant property of the current is governed by the area under the FS. In 
Fig. 6 the currents are shown for an elliptical FS with (3>«. Different insert Fig. 6 
current characteristic in this figure correspond to Kp increased up to 
2% from its initial value. The plot shows that in this way the paramag- 
netic currents can be converted into diamagnetic and vice versa. The 
characteristic is therefore an almost periodic function of the volume 
under the FS. This important property is reminiscent of the current be- 
havior in a single channel case, when the number of electrons in the 
system changes by two [0, ||. Here we show that this property of a 
single channel survives the transition to the 2D case. 

So far we considered a perfect cylinder. However it has been shown [[I], 
|5|, [T2[ that the phenomenon will survive modest scattering, both elastic and 
inelastic. Let us assume that our mesoscopic system contains a small number 
of impurities (ballistic regime). The average current, where the average is 
taken over impurity configurations, has been calculated in |[12|| . The formula 
for the total average current for 2D cylinder is 




(17) 



where A is the mean free path. For a sufficiently clean material A is of the 
order of a few microns and in the mesoscopic regime the impurities do not 
decreases the current significantly. 



4 Spontaneous Currents in the Mesoscopic 
Cylinder 
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field in the ID system of the mesoscopic size in T > is replaced by 

^'f^ ^ ^A)*^), (18) 
fri ttt* i _ exp v O y 

where T* is the characteristic temperature, defined by the energy gap between 
energy levels at the Fermi surface 

T =^- (19) 



The formula (|18|) is valid assuming the grand canonical ensemble [0] . It is ap- 
plicable in our case because single channels can exchange electrons according 
to K x , which determines their chemical potential. Temperature influences 
the current in such a way, that all discontinuities in the I— (j) characteristic 
(all transitions in the current value) are smoothed, and the maximum of the 
amplitude decreases. 

Consider now a 3D system of the cylinder geometry made of a clean non 
superconducting material. The circumference, height and the width of the 
wall are denoted by L x ,L y and L z respectively. We assume that the width 
of the cylinder is small compared to the other dimensions L z <C L x , L z <C L y 
and therefore we can assume with a good approximation that the vector 
potential A does not depend on z. 

Generalizing equation (||) to a 3D system the current can be obtained by 
replacing K x in (|1^) by 

K x = aK F " 
where a(3^ = 1 and 



1 (ky(r)Y (Us) 
\ \(3K F ) {jK F 



(20) 



S7T 

k z (s) = — } s = 0,1,2,.... (21) 



L 



z 



Therefore we obtain 



M P oo AT T 
r=0 s=0 1=0 " ± \ 



1 - 



'k y (r)V 



k z (s 



exp 



exp {-^] 



x 
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X COS 



lKpL x a r 



1 (kM 




sinf^V (22) 



Changing the Fermi surface in the same way as in 2D case yields a very similar 
effect. In cases where the current coherence was high the third dimension 
substantially increases the number of channels (P ~ L z ) and enhances the 
current. 

It has been recently shown in a series of papers |2|, |J that a metallic or 
a semiconducting system made of a set of mesoscopic quasi ID rings stacked 
along certain axis can exhibit a transition to a low temperature state with a 
spontaneous orbital current. 

In this chapter we concentrate on searching for a possibility of obtaining 
such spontaneous currents in a clean metallic or semiconducting 3D cylinder. 
We consider here the case when the magnetic field (f> enclosed by the cylinder, 
which drives the persistent current I ((f)) around it, is the sum of the externally 
applied flux (f> e and of the flux (f>i from the persistent current itself |2|] 

<P = (f>e + CI, (23) 

where C is the self- inductance of the cylinder (for cylinder with height L y , 
width L z and a radius R = f^, C = &^(^ L 2 y - R? - R)). Most of the- 
oretical discussions [7] neglect the second term in ([23]). Relations (53) 
and fl2"2"| ) represent the system of two selfconsistent equations to calculate the 
current I ((f>) and the magnetic flux (f>. The question of existence of the spon- 
taneous persistent currents reduces into a problem whether these equations 
have stable, non- vanishing solutions at e = 

2 = (24) 



If there is any nonzero solution of (|24[) then we have the spontaneous current, 
corresponding to the point of intersection the straight line (LHS of (p4|)) and 
the nonlinear curve (RHS of (|4j)) which is given by fl22|) . The nonzero 
solution is stable provided the intersection is at the place where the slope of 



the curve (|22[) is negative. 

We have made model calculations for a cylinder with the following pa- 
rameters L x = 25000vrA, L y = lOOOOA, L z = lOOA, E F = 7eV with a = (3 = 



11 



7 = 1 and n = 4, 5 (Figs 7, 8). We see that for n > 4 there is always a insert Figs 7, 

nonzero stable solution at sufficiently low temperatures (circles at Figs 7 and 8 

8). The temperature at which the spontaneous current occurs increases with 

increasing n. One can show also that the condition a — (3 — j — 1 is not 

necessary and the spontaneous currents occur even if the above parameters 

are different. It means that for all the Fermi surfaces similar to rectangular 

the spontaneous currents can be obtained. 

Similar transition to a spontaneous currents occur for a generalized tri- 
angular FS i.e. for n = 1. In this case the FS is of a tetrahedral shape and 
the geometrical condition necessary to obtain spontaneous currents is 

^aL x = fi(3L y = wyL z , (25) 

where /i and v are positive integers. 

For 2D systems we do not observe spontaneous currents neither for high 
value of n nor for reasonable values of a and (3 coefficients. 



5 Conclusions 

In this paper we presented same model considerations of persistent currents 
in clean multichannel systems of a cylindrical geometry (ballistic regime). 

We discussed persistent current behavior due to the pure Aharonov-Bohm 
effect — we neglected attenuation of the current due to the field penetrating 
through the system. We also ignored the spin-orbit coupling. One can es- 
timate that such assumptions @ pose no serious problems for the cylinders 
under study. We also neglected here the Coulomb interaction. Recently it 
has been shown that the Coulomb interaction does not influence the persis- 



tent current behavior in clean systems [13|, whereas it enhances the current 



in the diffusive regime |L4]]. There is also an indication that the Coulomb 
interaction of the Hubbard type enhances the typical value of the persistent 
current in disordered rings in any dimensions fll~5|. The decrease of currents 



with disorder was taken into account (see eq. |I7]) following [12 

The current in multichannel system depends strongly on the phase cor- 
relations between currents of different channels. We discussed persistent 
currents as a function of a shape of the Fermi surface by comparing them to 
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the standard circular FS. In case of the elliptical FS the current-flux charac- 
teristics show substantial increase of the current magnitude with increasing 
and attenuation with decreasing curvature of the Fermi surface at the point 
where it crosses the k x axis. The triangular FS enhances the current pro- 
vided its parameters obey the geometrical condition ([TBI) . The departure 
into concave FS only slightly increases the current. In the opposite case of 
the convex FS the current amplitude increases from its value in case of the 
circular FS to the maximal value in case of the rectangular FS (n — > oo). We 
found also that the current-flux characteristic is an almost periodic function 
of an area under the FS. 

We also investigated the possibility of spontaneous currents in 2D and 
3D cylinders. We found that for 3D cylinders spontaneous currents can be 
obtained for n > 4 and for n = 1 provided the geometrical dimensions of the 
cylinder are properly adjusted. This can be a hint for a possible experimental 
verification of the currents. The characteristic temperatures T c increase with 
increasing n > 4. The most favorable situation takes place for n — > oo 
what is equivalent to perfect correlation of currents from different channels. 
Such Fermi surfaces can be obtained e.g. for bcc crystals in tight binding 
approximation for nearly half filled or half filled band. Preferable materials 
would be metals or semiconductors with a long phase coherence length and 
with a small number of impurities. 

6 Acknowledgments 

Work was supported by Grant KBN PB 1108/P03/95/08. E.Z. is grateful to 
Prof. W.Brenig and Prof. W.Zwerger for useful discussions. 



13 



7 Figure captions 



Figure 1. Persistent current vs flux at T=0, for elliptical Fermi surfaces 
in = 2), with /3 = 0.1 (dashed line), (3=1 (solid line) and (3 = 10 (dotted 
line). In the inserted circle the magnified curve (a = (3 = 1) show the details 
of the current behavior. 

Figure 2. Convex Fermi surfaces according to Eq.(H) for a = (3 = 1 and 
n = 1, 2, 5 and oo. 

Figure 3. Concave Fermi surfaces according to Eq.([|) for a = (3 = 1 and 
n = 0.8,0.6 and 0.4. 

Figure 4. Persistent current vs flux at T=0 and a = (3 = 1, for spherical, 
n = 2 (solid line), triangular, n = 1 (dashed line) and concave n = 0.4 
(dotted line) Fermi surfaces. 

Figure 5. Persistent current vs flux at T=0 and a = j3 = 1, for convex Fermi 
surfaces, n = 2 (solid line), n = 5 (dotted line) and rectangular FS, n — > oo 
(dashed line). 

Figure 6. Persistent current vs flux at T=0 and (3 ^> a and n = 2, for the 
elliptical Fermi surface and different volumes Kp l < Kp 2 < . . . < Kp 5 . The 
characteristic is an almost periodic function of Kp (only half of the period 
in Kp is shown). 

Figure 7. The graphical solution of a set of selfconsistent equations (^2j) 
and (|24l) for different temperatures and the convex Fermi surface n = 4, 
a — (3 — 1. The nonzero crossings of the straight line (|24l) with the current- 
flux characteristic (|22D denoted by circles correspond to spontaneous currents 
in the 3D cylinder. 

Figure 8. The graphical solution of a set of selfconsistent equations 
and (|24l) for different temperatures and the convex Fermi surface n = 5, 



a — (3 — 1. The nonzero crossings of the straight line (|24|) with the current 



flux characteristic (p2|) denoted by circles correspond to spontaneous currents 
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in the 3D cylinder. 
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